Universal spin-polarization fluctuations in ID wires with magnetic impurities 
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We study conductance and spin-polarization fluctuations in one-dimensional wires with spin-5/2 
magnetic impurities (Mn). Our tight-binding Green function approach goes beyond mean field 
thus including s-d exchange-induced spin-flip scattering. In a certain parameter range, we find 
that spin flip suppresses conductance fluctuations while enhancing spin-polarization fluctuations. 
More importantly, spin-polarization fluctuations attain a universal value 1/3 for large enough spin 
flip strengths. This intrinsic spin-polarization fluctuation may pose a severe limiting factor to the 
realization of steady spin-polarized currents in Mn-based ID wires. 

PACS numbers: 72.25.-b, 73.63.-b, 73.63. Nm 
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I. INTRODUCTION 

Spin-related effects in novel solid state heterostruc- 
tures give rise to a rich variety of fascinating physi- 
cal phenomena. These spin-dependent properties also 
underlie a potential technological revolution in conven- 
tional electronics. El This new paradigm is termed "Spin- 
tronics" . A particularly interesting theme within this 
emerging field is spin-polarized transport in semiconduc- 
tor heterostructures. This topic has attracted much at- 
tention after the fundamental discovery of exceedingly 
long spin diffusion lengths in doped semiconductorscl fol- 
lowed by the seminal spin injection experiments in Mn- 
based heterojunctions.B 

Theoretically, a number of works have addressed is- 
sues connected with spin-polarized transport- These in- 
clude, for instance: spin filtering,Q spin waves,Qand quan- 
tum shot noiseja - all in ballistic semimagnetic tunnel 
junctions - and mesoscopic conductance fluctuations in 
Rashba wires. Spin-dependent phenomena in connection 
with localization effects should bring about exciting new 
physics. 

Here we investigate conductance and spin-polarization 
fluctuations for transport through one-dimensional wires 
with spin-5/2 magnetic impurities, e.g., Mn-based II- 
VI alloys such as ZnSe/ZnMnSe/ZnSe. The exper- 
imental feasibility of these wires has already been 
demonstrated!™ In these systems, the conduction elec- 
trons interact with the localized d electrons of the Man- 
ganese via the s-d exchange coupling.t2l UCF in Mn- 
based submicron wires was first experimentally studied 
in Ref. ||. We describe transport within the Lan- 
dauer formalismliil and calculate the relevant transmis- 
sion coefficients via non-interacting tight-binding Green 
functions £j 

We treat the s-d interaction beyond the usual mean- 



field theory thus accounting for spin flip scattering. In 
a certain parameter range we find that &nin-flip scat- 
tering suppresses conductance fluctuations!^ (below the 
UCF value for strictly ID wires) while enhancing the cor- 
responding spin-polarization fluctuations. More impor- 
tantly, we show that the spin-polarization fluctuations 
attain a universal value ((S() 2 ) = 1/3 for strong spin-flip 
scattering. This large spin-polarization fluctuation may 
pose a fundamental obstacle to attaining steady spin- 
polarized currents in Mn-based wires. 



II. MODEL HAMILTONIAN 

We consider a one-dimensional tight-binding chain, see 
Fig. |l|, of N spin s = 5/2 magnetic impurities coupled 
to ideal leads (sites n < 1 and n > N). We separate 
the electronic and impurity-spin degrees of freedom and 
treat the latter classically (static scatterers). The two- 
component electron wave function, vb — (ip^ipi) is then 
governed by the Schrodinger equation with a Hamilto- 
nian 
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Here, Ho is spin independent with clcmcntsEl 



2jO~nm — J$nm+1 — 7<5 nm _i + V n 5„ 
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where V n is the potential at site n and 7 = h 2 /2ma 2 , with 
a being the "lattice constant" . In the leads H itself gives 
rise to the usual dispersion relation e(k) — 27(1 — cos ka). 

In the following, a =\ = 1/2 and a =| = -1/2. We 
restrict ourselves to zero magnetic field so that the block- 
matrices H aa have elements given by 



{H aa } % 
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FIG. 1: One-dimensional tight-binding chain with N magnetic s = 5/2 impurities (mutually un-correlated, each spin equally 
distributed among the 6 spin states) coupled to ideal impurity- free leads (sites n < 1 and n > N). 



which is a Heisenberg-like interaction of the spin of the 
electron (a) with the z-component spin of the impurity 
S = (S x , S y , S z ). The off-diagonal block-matrix Hfi = 
if J I contains the interaction of the electron spin with the 
x and y-components of the impurity spins which leads to 
spin flip 



III. TRANSPORT PROPERTIES 



(4) 



We consider a sufficiently weak coupling between the 
impurity spins so that they can be considered mutu- 
ally uncorrelated, i.e., no magnetic ordering. The z- 
component of each spin is equally distributed among the 
6 spin states and the x and y-components are uniformly 
distributed with the constraint that S 2 = S 2 + S 2 + S 2 — 
s(s + 1), see Fig. [jj 

We study transport in the low-temperature linear re- 
sponse limit within the Landauer formalismEU 



(5) 



Here, T is a 2 x 2 matrix with the elements T aa i being 
the transmission probability of an electron from a state 
with spin a' in one lead to a state with spin a in the 
other lead. From Eq. (||) we now define the degree of 
spin polarization 



c _ h - h 



T TT + T T1 ~ T i1 
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which we will focus on in this paper. 

Green function method. The transmission matrix T is 
related to the retarded Green function 



G[e) = (e ■ 1 - H - -E(e))- 1 
via the Fisher-Lee relatiorB 



,{e) = [hv{e)f\{G^{e)} 



Nl 



(7) 



(8) 



where v = H 1 de/dk is the group velocity in the leads. 
In Eq. (0) the 2N x 2N matrix H is the Hamiltonian 
truncated to the N lattice sites with magnetic impurities. 
The effect of coupling to the leads is contained hi— the 
2N x 2N retarded self-energy matrix with elementalj 



{5w(e)}r 



-~fe* k ^ a 6 aa ,6 nm (6 ln + 5 Nn ). (9) 



N = 1 case. A chain with a single impurity is a simple 
illustrative example where analytical progress is possible. 
After performing the straightforward matrix inversion in 
Eq. (f7) we find 



C(e) 



-V\J Z S Z 



V 2 + e(4 7 - e) + (J • S) 2 



(10) 



where J — (J x , J y , J Z ) T ■ In zero magnetic field (S z ) = 
and (S 2 ) = 35/12. This implies that (C) = both with 
and without spin flip whereas the fluctuations are finite. 
The analytical averaging is of course complicated by the 
presence of S z in the denominator, but for isotropic cou- 
pling J x = J y = J z = J , we have (J • S) 2 = Jqs{s + 1) 
so that S z only shows up in the numerator, i.e. 



((Kf) 



35 



V 2 T 2 



12 [V 2 + e{A 1 - e) + J 2 s{s + l)f 



(11) 



In the absence of spin-flip (J x = J y = 0) the fluctua- 
tions are enhanced due to the replacement of s(s + 1) — > 
S z < s(s + 1) in the denominator (the final expression 
for the fluctuations is much more complicated) and this 
means that spin-flip will lower the fluctuations of £! Of 
course this trend is strictly valid for N = 1 , but in a lim- 
ited parameter range, this trend is still true for larger N 
values. 

Finite N case. For a finite number of impurities the 
problem is not analytically tractable and we study the 
problem numerically by generating a large ensemble (typ- 
ically 10 5 members) of spin configurations. For each spin 
configuration we calculate Eqs. (0,||) numerically. In our 
simulations we use the following parameters: Ef = 7, 
J z = 7/2, V n — (i.e., we neglect spatial disorder), and 
varying spin- flip coupling strengths ^ J x — J y ^ 7. 



IV. RESULTS AND DISCUSSIONS 

Figure || shows the distributions P(C), P(T aa -), and 
P(T aa i) for N — 10 and increasing strengths of the spin- 
flip coupling J x = J y . The distribution P(Q is symmet- 
ric around £ = which implies that on average there is 
no spin filtering, (£) = 0. The distribution P(C) first gets 
narrower for spin flip in the [0, 0.15 7 ] range (not shown) 
and then broadens as spin flip further increases. For suf- 
ficiently strong spin-flip scattering the distribution ap- 
proaches that of the uniform limit in which P(C) = 1/2. 
In this limit P{T aa ) and P{T aa i) coincide and so do all 
average transmission probabilities (T aa i). As we discuss 
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FIG. 2: Distributions P(C), P(T„a), and P(T CTCT <) for different 
spin-flip scattering strengths = J M in the case of N = 10. 
The dash-dotted line in the lowest panel indicates the uniform 
limit P(C) — 1/2 [note the magnification of -P(C) by of factor 
of 5] attained for strong enough spin-flip scattering. 



below, the initial narrowing and subsequent broadening 
of P(C) with spin flip gives rise to a minimum in the 
fluctuation of £, Fig. ||. 

Universal spin-polarization fluctuations. In the limit 
of a short spin-flip length <„ « L wc in general find a 
uniform distribution -P(C) = 1/2, Fig. 0. This uniform 
distribution yields the universal value ((<5£) 2 ) = 1/3 for 
the spin-polarization fluctuations. Figure || clearly shows 
that this universal value is attained for increasing spin flip 
strengths and is indeed independent of N. Interestingly, 
Fig. Halso shows a minimum at around J x = J y = O.I57. 
This minimum can be attributed to two competing en- 
ergy scales: the longitudinal (~ J z ) and the transverse 
(~ J x , J y ) parts of the s-d exchange interaction, Eqs. (||) 
and (JI]), respectively. A simple "back-of-the-envelope" 
calculation shows that these two competing scales are 
equal for J x = J y = J z /J2s(s + l)/3j = O.2O87. The 
vertical dashed line in Fig. [^indicates this value. Observe 
that (((5£) 2 ) becomes larger for increasing N. This hap- 
pens because P(C) broadens for larger TV's (the travers- 
ing electrons see a wider region with random spins). This 
is similar to the broadening due to increasing spin flip 
strength. 



We should mention that the distribution P((), and 
consequently ((8() 2 ), change dramatically for Ep < J z . 
In this regime, P(Q becomes [/-shaped (not shown) be- 
cause of the dominant filtering due to the "end states" 
with Sj tZ — ±5/2. This qualitatively different P(£) yields 
a monotonically decreasing ((^C) 2 ) as a function of spin- 
flip strength. Here the universal ((<5£) 2 ) = 1/3 value 
is approached from above for large spin flip strengths 

(Jx = Jy~ l)- 
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FIG. 3: Average fluctuations ((^C) 2 ) as a function of spin- 
flip strength for N = 10, 20, and 30. The dashed horizon- 
tal line indicates the universal value 1 / v3 obtained from the 
uniform limit P(C) = 1/2. The vertical dashed line indicates 
where the spin-flip rate is comparable to \J z a\/h. 



Suppression of conductance fluctuations. Whereas the 
fluctuations in the spin polarization £ remain finite in the 
strong spin-flip scattering regime (Fig. ||), we find that 
the fluctuations of the conductance g are strongly sup- 
pressed in this limit. This is illustrated in Fig. which 
shows the average conductance and its fluctuations as 
a function of spin- flip scattering for N = 10, 20, and 
30. Note that ({Sg) 2 ) 1 ^ 2 is much more sensitive to spin 
flip than (g). In addition, for all N we essentially have 
{(5g?) 1/2 > (9) for J x = J y -> and ((Sg) 2 ) 1 / 2 < (g) 
for J x — J y — > 7. Figure [| clearly shows the conductance 
fluctuations get suppressed for increasing N. The hori- 
zontal dashed line shows the UCF value (0.73/2 = 0.365, 
see e.g. Ref. [TH]) for a ID wire in the metallic regime. The 
spin-related conductance fluctuations do not approach a 
finite value for increasing spin flip scattering. It actually 
seems to go to zero. This is in contrast to the spin- 
polarization fluctuations, Fig. ||, which attain a universal 
value ((^) 2 ) 1/2 = 1/V3 for strong spin-flip scattering. 
Incidentally, we observe that {(Sg) 2 ) 1 / 2 and ((^C) 2 ) 1 ^ 2 
also present contrasting behavior for increasing N (and 
Ep > J z ): the former gets suppressed while the latter 
gets enhanced (cf. Figs. § and [I). 

Spin disorder as spatial disorder. To some extent, 
the s-d site interaction considered here plays the role 
of spatial disorder in the system with a mean-free path 
£j. Let us consider first the case with no spin flip (i.e., 
Jx = Jy = 0). In this case, the term J z o~S n ^ z acts as a 
"random" spin-dependent potential along the chain (here 
the site potential has some internal structure). As shown 
in Fig. |j the conductance fluctuations for zero spin-flip 
scattering are larger than, slightly above, and slightly be- 
low, the UCF value for TV = 10, N = 20, and N = 30, 
respectively. For increasing N we go from the metal- 
lic regime (L = Na <C £j z ) with vanishing fluctuations 
and a Gaussian P(g) strongly peaked near g ~ 2e 2 /h 
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FIG. 4: Average conductance (<?} and its fluctuations 
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as a function of the spin- flip scattering strength for 



N = 10, 20, and 30. The conductance fluctuations are much 
more sensitive to spin flip than the average conductance: the 
former gets strongly suppressed for increasing spin flip rates. 



to the strongly localized regime (L £j r ) where it is 
well-known that P(g) is strongly peaked near g ~ 
with a log-normal distribution so that fluctuations can 
be comparable to the mean value.ll3 This is in accor- 
dance with numerical studies with different continuous 
distributions of the "on-site" potential (e.g. Gaussian 
or uniform distributions). E3 In Fig. 4 the "small" mean 
values (g), for N = 10,20, and 30, indicate the onset 
of localization with fluctuations comparable to the mean 
value. As N gets larger-conductance fluctuations are as 
expected suppressecOia. 

Role of spin-flip scattering. Spin flip clearly suppresses 
conductance fluctuations, Fig. [|. This can be under- 
stood from Eq. (4) being a complex number with a ran- 
dom phase which makes spin flip act as a source of "de- 
coherence" [the total wave function is, of course, fully 
coherent]. Furthermore, spin flip mixes all the S n _ z com- 
ponents on each site thus smoothing the potential seen by 
the traversing electron and hence reducing conductance 
fluctuations. This is true for both ef > J z [except for the 
window [0,0.157] in which P(Q narrows] and Ef < Jz- 



"Truly" universal fluctuations. Why is ((5Q) 2 ) 1 / 2 uni- 
versal even for short spin-flip lengths £ a <C L (strong 
spin-flip scattering) while ((Sg) 2 ) 1 / 2 is clearly suppressed 
below the usual UCF value in this limit? It is well 
known that conductance fluctuations are suppressed in 
the incoherent limit More specifically, in ID wires with 
itp <C L, £ v is some "dephasing length", the suppression 
factor is y / L/£ v (see Ref. Q). Interestingly, we can 
likewise understand the suppression of ((Sg) 2 ) 1 / 2 seen in 
our simulations by viewing spin-flip .scattering as pro- 
ducing "dephasing" with £ v ~ £j x tB. For the spin- 
polarization fluctuations, however, the picture is slightly 
different: here we divide our system in Nl — L/l v seg- 
ments. To each of these we can associate an average spin 
polarization (Q) = (i: 1..Nl) and a corresponding 
spin-polarization fluctuation ((SQ) 2 ). Neither (Q) nor 
((S(i) 2 ) are additive quantities like (g) and ((Sg) 2 ) ("ex- 
tensive versus intensive" properties). Sensible global av- 
erages for the whole system are then £ = J2i(d) — 

and (S() 2 = T^J2 t ((SQ) 2 )- We should expect (S(f = 
((S(i) 2 ) = ((S() 2 ) if the- system is ergodic. Hence uni- 
versal spin-polarization fluctuations are not suppressed 
for large spin-flip scattering in contrast to conductance 
fluctuations. 



V. CONCLUDING REMARKS 

Spin-flip scattering in Mn-based wires reduces conduc- 
tance fluctuations while enhancing spin-polarization fluc- 
tuations in a limited parameter range. Remarkably, spin- 
polarization fluctuations reach a universal value 1/3 for 
large spin-flip scattering in which the conductance fluctu- 
ations vanish. This universal value should manifest itself 
in time- and polarization-resolved photo-luminescence 
measurements. More important, these sizable spin fluc- 
tuations may limit the possibilities for steady spin injec- 
tion in these systems. 
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